A consistent quantum model for continuous photodetection processes 



o 
o 



> 
On 
00 
O 
t> 

o 
m 
O 

9 L,: 
^— > ■ 

G ■ 

s : 
cr 



M C de Oliveira*, S S Mizrahi^ and V V Dodonov* § 
Departamento de Fisica, CCET, Universidade Federal de Sao Carlos, 
Via Washington Luiz km 235, 13565-905, Sao Carlos, SP, Brazil. 
(February 1, 2008) 

We are modifying some aspects of the continuous photodetection theory, proposed by Srinivas 
and Davies [Optica Acta 1981 28 981], which describes the non-unitary evolution of a quantum 
field state subjected to a continuous photocount measurement. In order to remedy inconsistencies 
that appear in their approach, we redefine the 'annihilation' and 'creation' operators that enter 
in the photocount superoperators. We show that this new approach not only still satisfies all the 
requirements for a consistent photocount theory according to Srinivas and Davies precepts, but also 
avoids some weird result appearing when previous definitions are used, 
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I. INTRODUCTION 

The subject of quantum measurements is as old as the 
very foundation of quantum mechanics. For a long time 
the scheme proposed by von Neumann [1] has been preva- 
lent. According to this scheme, there are two kinds of 
evolution of quantum states: a unitary evolution obeying 
the rules of quantum dynamics in the absence of measure- 
ments, and an 'extra-quantum' dynamics resulting from 
a measurement, when the state vector suffers a sudden, 
instantaneous and irreversible transformation (or reduc- 
tion, or collapse) to one of possible eigenstates compati- 
ble with the observable measured. 

When one considers, for instance, a photocount pro- 
cess, the radiation impinges a photomultiplier tube and 
each read burst of electrons (current) is viewed as a man- 
ifestation of a single photon. A sequence of bursts is asso- 
ciated to photocounts. A classical theory describing this 
process was proposed by Mandel [2], however it resulted, 
under certain particular circumstances, in negative prob- 
abilities of counts. Thence refined quantum photocount 
theories were developed by Mandel, Wolf and Sudarshan 
[3,4], Glauber [5], Kelley and Kleiner [6], Mollow [7], 
Scully and Lamb [8], and others (see, e.g., the review 
[9] for more references). 

However these theories relied on the assumption of in- 
stantaneous measurement, although in reality photons 
are counted sequentially, one by one. In the 1980's Srini- 
vas and Davies [10] developed a theory describing actual 
photocounting events, which was adopted for estimat- 
ing the outcome in several problems, such as quantum 
non-demolition measurements [11,12], determination of 
field states under continuous photodetection process [13], 



quantum theory of field-quadrature measurements [14], 
conditional generation of special states [15], and for the 
control of the amount of entanglement between two fields 
[16]. 

Srinivas and Davies (SD) theory considers photodetec- 
tion as a continuous measurement, with no reference to 
a 'meter state'. The main quantities to be calculated in 
this theory are probability distributions for counts (or 
no counts). It is based on the assumption that in an 
infinitesimal time interval r only two processes may oc- 
cur: either one-count, characterized by a superoperator J 
acting on the field density operator p, or no-count, char- 
acterized by another superoperator, S T . Superoperator 
J in the SD theory has the form 



J p = ^apo) 



(1) 



where 7 is the detector efficiency and a (a)) is the field 
'annihilation' ('creation') operator ([a,at] — 1). Such a 
choice is based on the assumption that operator a sub- 
tracts one photon from the field. It is assumed [13] that 
just after one count (conceived as subtraction of a single 
photon from the field) the system state is given by 



p(t 



Jp(t) 
Tr[Jp(t)Y 



(2) 



where t + stands for t plus as infinitesimal time after. 
(From another point of view, the states described by 
means of the statistical operators of the form (2) were 
considered in Refs. [17-19] under the name 'photon- 
subtracted states'.) 

However, one can easily check that the mean number of 
photons in the state (2) (i.e., after counting one photon) 
equals 
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n(t + ) =n(t) + 



An 2 (t)-n(t) 
W) 



n(t) + q, (3) 



where An 2 



n 2 and q is Mandel's q parameter 



[20] characterizing the type of photon statistics in the 
initial state of field: for q < (q > 0) the field statis- 
tics is said to be sub-Poissonian (super-Poissonian) and 
for q = the statistics is Poissonian. However, equation 
(3) clearly shows that the mean number of photons may 
increase after one count, if the field statistics is super- 
Poissonian. Otherwise, the mean number will remain the 
same if the statistics is Poissonian, and will decrease only 
if the statistics is sub-Poissonian. For example, for a field 
represented by a Fock state |ra), one gets n(t + ) = m — 1 , 
exactly one photon less, whereas for any other field state 
this is not true. This point received special attention in 
[21]. Thus, we see that in general the common choice for 
J docs not really correspond one count to one less pho- 
ton in the field. Similar observations were made, e.g., in 
[13,18,22,23], however, without attempts to modify the 
theory of photocounting processes. 

Besides, Srinivas and Davies themselves [10] perceived 
that the superoperator J does not satisfy assumption 
(V) of their theory, namely the boundedness property 
Tr(Jp) < oo. In fact, J is an unbounded linear transfor- 
mation and consequently the counting rate is unbounded, 
thus not defined for all possible states. This fact leads to 
an ill-defined coincidence probability density. We shall 
return to this point in the next section. 

In this paper we propose some modifications in the SD 
photocount theory in order to satisfy all the assumptions 
proposed by its authors. Our motivation finds ground on 
the recent discussion about the role of the 'annihilation' 
operator a in quantum optics [21], since state a\ip) is not 
always a state whose mean number of photons is less than 
in Depending on the field statistics, state a\ip) may 
show much higher mean number of photons than state 
\tp). In [21] the authors suggested that instead of a and 
a' the exponential phase operators E- and E + should 
be considered as real 'annihilation' and 'creation' oper- 
ators in the photocounting theory. The introduction of 
these operators in the continuous photocounting theory, 
besides eliminating inconsistencies in the SD proposal, 
leads to new interesting results related to the counting 
statistics. 

This paper is organized as follows. In section 2 we 
briefly revise the main aspects of quantum photodetec- 
tion theories. In section 3 we present our model based 
on exponential phase operators. In section 4 we consider 
the field state evolution under continuous monitoring, 
but when no information about the number of counted 
photons is read out, or the 'pre-selection' state evolu- 
tion. Also, in this section we solve the master equation 
generated by the exponential phase operators instead of 
the annihilation/creation ones (at zero temperature) and 
consider several important special cases. In section 5 we 
discuss a physical meaning of the results obtained, point- 



ing at the principal differences in the predictions of the 
SD theory and our model, which could be verified exper- 
imentally. 



II. FUNDAMENTALS OF CONVENTIONAL 
QUANTUM PHOTODETECTION THEORIES 

The first quantum photodetection theory, developed 
independently by Mandel et al [3,4], Glauber [5], and 
Kelley and Kleiner [6] as a simple extension of the classi- 
cal theory [2], gave the following probability of k counts: 

P(M)=Tr{p:i[ 7 i/(*)] fe e-^«: 

= E(*)(l-7tr*(7«) fc ft.. (4) 

n=k ^ ' 

Here p is the statistical operator of the field (we consider 
a simplified model of a one- mode field), 7 is the detec- 
tor efficiency, I(t) is the average field intensity, :: stands 
for operator normal ordering and p n — (n|p|n) is the 
probability to have n photons in the given field mode. 
However formula (4) becomes obviously meaningless if 
jt > 1, when it can result in negative probabilities or 
unlimited mean number of counted photons as t — > 00. 
These troubles were removed in studies [7,8,10], whose 
authors, using different approaches, arrived at the same 
result, which consists, from the formal point of view, in 
the substitution yt — > 1 — exp(— -ft) in equation (4): 



P SD {k,t) 



OO / 

n—k 



(l-e-^) k (e-^r- k (n\p\n). (5) 



The SD photocount theory does not refer to a specific 
detector state. It was built by considering two kinds of 
events represented by superoperators acting continuously 
on the field state. The first one, represented by J, is a 
single instantaneous count event, while the other, repre- 
sented by St, is a no-count event for a time interval t. 
Thus, the operation (a superoperator) 

N t (k) = f dt k f dtk-i ■ ■ ■ 



Jo 

*2 



t>Z2 

x / dtiS t - tk JS tk -t k _ 1 ■ ■ ■ JS tl 
Jo 



(6) 



stands for the count of k photons from the field. The 
operator (6) projects continuously the initial field state 
p, and P(k,t) = Tr[N t (k)p] is the probability of counting 
exactly k photons in a time interval t. 

The assumptions or properties of the SD theory are as 
follows. 

(I) p — > N t (k)p is a linear positive map on the space 
T(7i) of trace class operators on the Hilbert space 7i, 
such that for any positively definite p satisfying Tr p = 1 , 
one has 
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< Tr[N t (k)p\ < 1. (7) 

(II) The sum over all possible counted photons satisfies 
the normalization condition Tr[T t/ o] = 1 where 

T t = Y J N t (k) = S t + T t _ t ,JS v dt'. (8) 

(III) Semigroup associative property 

N tl+t2 (k)= N tAk2)N tl {ki). (9) 

k\+k 2 =k 

(IV) Identity 

limJV t (0)p = p. (10) 

(V) Assumption of bounded interaction rate: there exists 
a number K < oo such that 

oc 

5^Tr[JV t (fc)p] < Kt, for i > 0. (11) 

k=l 

(VI) Assumption of ideality: the operation 

S t = N t (0) (12) 

transforms pure states into pure states, noting that Stp 
is a pure state if p is a pure state. 

As a matter of fact the process of photocounting by 
a macroscopic detector is the resultant of many micro- 
scopic fundamental interactions between the electromag- 
netic field and atoms composing the detector. But how 
these fundamental interactions can be specified by a rea- 
sonable model for the count (or no-count) jump super- 
operator in a phcnomenological macroscopic photodetec- 
tion theory? It is worth citing the corresponding words 
from the Srinivas and Davies's paper [10]: " As a simple 
model for the measurement performed by the photode- 
tector, where an n-photon state will be converted to an 
(n — 1) -photon state whenever a photon is detected, it is 
reasonable to set" the single photon counting superoper- 
ator J in the form (1). We see that this operator was not 
derived from fundamental processes, but set ad hoc due 
to its apparent simplicity and reasonability. 

As soon as J is chosen, the superoperator for the state 
evolution between consecutive counts can be derived as 

S tP = e Yt pe YU , Y = -iH - (13) 

once it conserves the probability in a regular point pro- 
cess: 

Tr[Jp + Yp + pVt] = 0. (14) 

Note, however, that the superoperators J and S t de- 
fined above do not satisfy condition (V) , since 

Tr(Jp) = 7 Tr(pat a ), (15) 



or more generally, 

Tr[J k p\ = 7 fc Tr[ /9 : n k :] (16) 

is unbounded, so, not defined for all states. The viola- 
tion of assumption (V) prevents a consistent definition 
for all states of some important functions. For example, 
although one can define the elementary probability den- 
sity of counting k photons in a time interval t, attempts 
to calculate the coincidence probability density 

h(t 1 ,...,t k ) = Tr{T t _ tk J...JT tl p) (17) 

of counts observed at each of the times t\,...,tk together 
with other possible counts in between, when the detector 
is making measurements for a time interval t, result in 
serious problems in the SD theory. 

One of the goals of our paper is to show how one can 
avoid the violation of assumption (V) using other opera- 
tors instead of a and a 1 *. 



III. QUANTUM COUNTING PROCESSES WITH 
EXPONENTIAL PHASE OPERATORS 

Our idea is to use, instead of operators a and a 1 *, the 
so-called exponential phase operators 

E_ = (aWl)- 1 /^, (18) 
E+ ^aVa + l)- 1 / 2 , (19) 

introduced, as a matter of fact, by F London at the dawn 
of quantum mechanics [24] , although their systematic use 
began only after the paper by Susskind and Glogower [25] 
(for the history see [26] ) . The commutator of E_ and E + 
is the vacuum state projector, 

[E-,E+] = |0> <0| = A . (20) 

The normal ordered product of these operators is a com- 
plementary projector 

E+E- = 1- A = A = A f , A 2 =A. (21) 

Applying _E_ and E + on the number states one gets 

E-\n) = \n-\), E+\n) = \n+l) (22) 

with £L|0) = A|0) = 0, and [A, n] = (where n = a^a is 
the number operator). A useful property is 

e aA = A + e°A. (23) 

For other properties and generalizations see, e.g., [27-36]. 
Finally, it is worth to recall that E_ has as eigenstate the 
'coherent phase state' 

oo 

|^) = v ^ _ ^F^^|n), \z\ <1 (24) 

n=0 

as introduced in [37] and studied in [38-47] (see also [48] ) . 
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A. One-count event 



We redefine the one-count operator equation (1) as 



Jp = E-pE^ 



(25) 



Now J is bounded operator and the system state imme- 
diately after the 1-count process in the time interval [0, t) 
is transformed into 



Jp(t) _ Jp(t) 



Tr[Jp(t)] 1-po 



(26) 



where p = (0|p(i)|0) is the probability for the vacuum 
state. (Note that pure states E™\ip) were considered 
in another context in [49]. Mixed shifted thermal states 

Pth 11 ^ = E+PthE™ were studied in [50], whereas meth- 
ods of generating such states in a micromaser were dis- 
cussed in [51].) 

The mean number of photons in the state p(t) (26) is 



»(t + ) = J® 



Po 



1, 



(27) 



so, whenever a state p has none, or very small, contri- 
bution from the vacuum state, the counting operation 
extracts exactly one photon from the system, indepen- 
dently of the field statistics. For example, for the num- 
ber state p = |m)(m| (m ^ 0), n = m — 1, and for the 
coherent state p = \a)(a\ (a / 0), h = n/ (l — e~ n ) — 1, 
with n = \a\ 2 . 

On the other hand, for the thermal state 



1 



1 + n 



E 

n=0 



\ n 



l + n 



\n)(n\ 



(28) 



we obtain h = n, i.e., the mean number of photons is not 
changed, as expected. This is a correct description of 
a thermal system since taking out a single photon from 
a reservoir should not change its average number. Note 
that using the SD definition, Jp = apa) , one obtains the 
weird result n = 2n. So one perceives that using a and a ' 
for constructing a continuous photocount measurement 
leads to some inconsistent results. 



B. No-count event 

The time evolution between sequential counts is repre- 
sented by St = N t (0), a superoperator defined in terms 
of ordinary Hilbert space operators 



S t p 



Yt Y^t 

e Yt pe Y *, 



(29) 



where B t — e Yt is a semigroup element given in terms of 
the generator Y. The deduction of S t is conditioned to 
the relation 



Tr[Jp] = Tr[pR], 



(30) 



where R is the rate operator, related to Y by Tr(pR) = 
Tr(Fp + pY^), which substituted in (30) gives equa- 
tion (14). The theory requires that in the absence of 
counts the system has a unitary evolution, whose dynam- 
ics is governed by the free-field Hamiltonian H = fuva^a. 
Thus, the convenient choice satisfying (14) is 



Y = -iH-R/2 



-iH - -E+E. 
2 + 



(31) 



Taking into account equations (21) and (23), as well as 
the commutativity of operators H and A, one can easily 
calculate the result of action of nonunitary operator e Yr 
on a pure state (here r is an interval of time between 
counts) 



l^(r)> 



|^)H0hA)|0) + e -^/ 2 AhMT)), (32) 



where \4>h(t)} = exp(—iHr) is the freely evolved state 
vector (in the absence of measurements) . Thus the prob- 
ability of no- count event equals 

V (T) = \\\MT))\\=\W)\ 2 + e-^(mi>)- (33) 

For a mixed state, the same probability is given by 
Tr[SVp], and simple calculations result in the formula 



V (t) 



-Po(l 



"), 



(34) 



which, of course, coincides with (33) in the case of pure 
state. Note that lim r _ >00 'Po(7") = po, which means that 
the probability of no counts registered during an infinite 
time interval is equal to the probability of finding the vac- 
uum state in the measured state p. Formula (34) should 
be compared with analogous formula of the SD theory 
based on the operator Y of the form (13) 



Voir) 



n=0 



(35) 



Although equations (34) and (35) give the same limits for 
t — > oo, the intermediate time dependencies arc different. 



C. Continuous counting 

The continuous counting of fc-photons from a field in a 
time interval t is represented by a linear operator N t (k) 
acting on the system state during a time interval [0, t), 



P {k \t) 



N t (k)p(0) 
Tr[N t {k)p(0)Y 



(36) 



where p(0), or simply p, is the state of field prior to the 
counting process and P(k,t) — Tr(N t (k)p) is the prob- 
ability of counting k photons in t. The linear operator 
N t (k) can be written in terms of the operators St and J 
as 
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rt rt k rt 2 

N t (k) = dt k dt k -i--- / dh 
Jo Jo Jo 

xS t -t k JS tk -t k _ 1 ■ ■ ■ JS tl ■ 
Noticing however that 

JS t p = e-^U t (Jp), 

one gets 



St^JSt^t^ ■ ■ ■ JS tl = e-^U tk S t - tk J k 



where 



U t p = e- lHt pe lHt . 



(37) 
(38) 
(39) 
(40) 



It is convenient to introduce short notation for two 
partial sums of probabilities: 

k oo 

■Ah = ^2 Pn > Zk + 1 = E p n =l- A k . (41) 



n=0 



n=k+l 



Using (39) we can calculate the elementary probability 
distribution (EPD) of counts at the instants t\, t 2 , ■ ■ ■ , tk, 
if the total measurement time is t, 

V(t 1 ,t 2 ,-- -,t k ;t) = Tr [S t ~ tk JS tk - tk _^,- ■ ■, JS tl p] 

= 7 fe (e-^ k Pk + e~ 7 %+i) , (42) 

which should be compared to the EPD in SD theory 

V S D(ti,t 2 , ■■■,t m ;t) = 7 ™ m ! e -7(ti+t2+-+t m -mt) 



E 



17) 



e~ int Vr, 



In particular, for t = oo we obtain 

7>(ti,t2,-",*fc;oo) -7 fc e- 7tfc p fe , 



(43) 



(44) 



whereas the SD theory yields essentially different result 
for k > 2 

V SD (t u • • • , t m ; oo) = 1 m m\e-^+^+-+ t ^ Pm . (45) 

Only for the one-photon event, k = 1, both models pre- 
dict the same exponential probability distribution 

P(h) =7e- 7tl Pi =P(h)pi. 

We see that (42) corresponds to a Markovian process, in 
the sense that 

V(t 1 ,t 2 , ...,t k ; oo)/p fe = V(t x )V(t 2 - ti) ■ ■ ■ V(t k - t k -i), 

i.e., the EPD depends only on the last count, at time 
tk- On the contrary, in the SD theory the EPD depends 
on all times at which counts occur, moreover, each new 
count enters with increasing weight: 

V SD (ti,t 2 , - ■ ■ ,t m ;^)/p k =V{h)-2V(t 2 )- ...-mVitm). 



Nonetheless, both distributions, (44) and (45), have the 
same normalization 



ft k rt2 

dt k I dtk-i--- dt 1 V(t 1 ,t 2 ,- • • ,tk\°o) =Pk- 



o 



o 



We can write equation (37) as (k = 1, 2, . . .) 



(k)p = U t [ dt' 
Jo 



-it' 



it') 



k-l 



■ exp 



(fc-1)! 
x (J k p) exp [- |(t - i')A' 



(46) 



where t' = tk and the probability of occurrence of k 
counts in a time interval t is given by 



P{k 1 t)= Pk 



l-e" 7t E 



+ 



k\ 



OO 

Ep". 



(47) 



n—k 



which verifies the normalization condition 



£p(M) = i- 



(48) 



k=0 



The limiting value of equation (47) 



lim P(k,t) = (k\p\k)=p k 

t—>OQ 



(49) 



means that, asymptotically, the counting statistics coin- 
cides with the photon statistics (as it should), since we 
do not consider here the possibility of photons lost to the 
surroundings or the failure to count any photon exiting 
the cavity. If a photon leaves the cavity it is detected 
and counted, for sure. 

The moments of distribution (47) are given by 



k l t = J2P(k,t)k l = {k l ) 



k=0 



3=0 



3- 



k=j 



k=0 



k l (~ft) h 
k\ 



1=0 



where (k l ) = Y^=o ^Vk- In particular, for the mean 
number of photons one verifies that lim^oo k t = (k). 
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D. Examples of field states 

Now let us consider specific field states, comparing the 
probabilities of fc-counts resulting from the Srinivas and 
Davies formula (5) and from our formula (47). For a 
field initially prepared in a number state \m), m > k, the 
probabilities are 



ni 



P* u (k, t) = (£)(!- e -T*) fe (e- 7 *) m - fe , (51) 



P(k,t) = e 



-it 



k\ 



(52) 



Different probabilities of a fc-count reveal different 
physical schemes of counting. In SD's theory, P SD (k,t) 
is a binomial distribution reflecting an underlying one- 
dimension 'random walk' process: the factor (1 — e~ 7 *) fc 
is the probability of k photons leave the cavity, while 
(e _7 *) m k is the probability that photons stay in the cav- 
ity, thus being not counted. On the other hand, the dis- 
tribution (52) is Poissonian. That means that each pho- 
ton leaving the cavity is counted, resulting in the same 
counting statistics as for falling raindrops in a small area, 
roughly the drop size. The last process is the only one 
consistent with the initial assumption that every photon 
leaving the cavity is counted. In figure 1 we compare 
both probability distributions, (51) and (52), for a field 
with m = 5, and selected /c-counted photon numbers. 
P(k,t) for (52) reveals a more spread shape then (51), a 
characteristic of Poissonian processes. In that figure we 
included the probability distributions for k = m. The 
SD theory gives the expression P SD (m, t) = (1 — e~ 7 *) m , 
while the present model gives P(m, t) = <f> m _i(7i), where 



n=0 



n=k+l 



x'- 



(53) 



Notice that for k > m, both theories give a zero- valued 
probability distribution, which is a signature of the fixed 
number state for the field inside the cavity, meaning that, 
in any time interval t, it is not possible to count more 
photons than those present in the field at time t = 0. 
For the coherent state \a) SD theory gives 



P SD (k,t) = I [H 2 (l- e -r*)]*exp [-^l 2 (l-e" 7 ')] , 



while in our present approach we get 
P{k,t) = 1 [( 7 i) fe e- 7t $ fe _i(|a| 2 ) + H 2fc e-H 2 $ fc ( 7 t) 

In figure 2 we display, for comparison, both distribu- 
tions for a coherent state with average photon number 
| a | 2 = 5, and selected counted photon numbers. 

For the thermal state (28) we obtain the distributions 

[n (1 — ■• 



P SD (k,t) 



[l + n(l- e" 7 *)] 



fe+i ■ 



P{k,t)= e —(^ ' 



M \n+l) + (l + n)' 



,fe+i 



which arc displayed in figure 3, for n = 5. For both, 
the coherent and thermal states, the probability distri- 
bution for the present model also show a more spread 
shape than the distribution for the SD theory. Perhaps, 
this feature may be important in the distinction of both 
models for photocounting, by experimental evidence. By 
repeated experiments of photocounting the output of a 
leaking cavity, one may reconstruct those distributions, 
assuming that any other incoherent (dissipative) process 
is absent or negligible for the whole process. 



IV. 'PRE-SELECTION' STATE EVOLUTION 



Having no knowledge about the number of counted 
photons after a time interval t, the field state being con- 
tinuously monitored, but with no readout is given by 



~p(t)=T tP = Y^N t {k)p. 



(54) 



fc=0 



It is referred as the preselection state [52]. Summing 
over k in equation (46) we obtain 



p(t) = S tP + J2 N t(k) P - U t {e-^pe-^ 2 

'e Jt ' (Jp)] e-i A W 2 dt j . (55) 



fe=i 



e -7*' e -7A(t-t')/2 



The probability of having n ^ photons in the cavity 
after continuous measurement for time t is equal to 

p n (t) = (n\p(t)\n) - e- 7 ' £ ij fpn+i , (56) 



where p n +; is the probability at t = 0. The probability 
to have no photons, in the cavity, at time t is obviously 

OO OO / ,w 

po(t) = 1 = VU-A, (57) 



n=l 



1=0 



where Ai is defined in (41). So, if one waits very long 
time the cavity will eventually end in the vacuum state, 
all photons being absorbed by the detector and counted. 
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A. The master equation 

Another way to obtain formula (56) is to use the 
phenomcnological master equation of the Lindblad form 
(hereafter we suppress the tilde over the operator p) 



E + E_p(t)- p{t)E + E_). (58) 



It should be compared with the 'standard master equa- 
tion' [53,54] for the amplitude damping model used in 
the SD theory, derivable from the interaction of a single 
electromagnetic (EM) mode (or 1-D harmonic oscillator) 
with an environment made of many harmonic oscillators 
at T= OK (see, e.g., [55] for the most recent applications) 

9 -^ = l [2ap(t)<f - a)ap{t) - p(t)a*a] . (59) 

If the EM mode is within a dissipative cavity, equation 
(59) can be viewed as describing the field state of an 
uncontrollable "photon-leaking" process to the environ- 
ment, although it could be not necessarily one- by-one, 
since J = japa^ does not produce a minus-one-photon 
state. In contradistinction, if there is no vacuum, for 
sure, inside the cavity, equation (58) describes the pro- 
cess of subtracting photons from the cavity, sequentially 
and one-by-one. 

Equation (59) results in the equation for the mean pho- 
ton number 



number due to the nonconstant coupling parameter T, 
following the choice of operators for picking out a 
single photon from the field. 

Note that equations (60) and (62) practically coincide 
in the case of low field intensities, when p n <C 1 for n > 2 
(this is especially clear from the expressions in the form 
of series over p n ). On the other hand, the 'saturation' of 
the decay rate in the case of high field intensities can be 
also understood, if one takes into account a possibility of 
large dead time of the detector: in such a case, for short 
intervals of time, the detector can count (with some ef- 
ficiency) only one photon, independently of the number 
of photons in the cavity or beam. This example shows 
that the decay rate (62) by no means can be considered 
as 'unphysicaP apriori. 

It is immediate to see that the solution to equation 
(62) is 



oo oo 



(ata>(*) = ]T np n (t) = e~i* £n£ { ^fp n+l 

71=1 1=0 



71=1 



= e-' t 



.71=1 



= 1 71=1 
OO 

+ (7*) H — 



71=1 



(64) 
(65) 



-d(a) a) (o) / 'dt = 7(a t a) = 7^np„, 



71=1 



whose solution 



(60) 



(61) 



shows that the rate of decrease of the mean number of 
photons is proportional to the mean number of photons 
present inside the cavity (or in the beam, whether beams 
are considered) resulting in an exponential decay as time 
goes on. 

Equation (58) leads to a quite different differential 
equation for the mean photon number 



-d(ata)^V^ = 7(l-f>o)=7$>" 



(62) 



71=1 



according to which, the rate of change in the mean num- 
ber of photons is proportional to the probability that 
there are photons in the cavity, independently of their 
mean number. If initially po = 0, this means that there 
are photons in the cavity for sure and equation (62) be- 
comes 

-d(a^a)^/dt\ t=o = r (Va>|(f>) ( a ta>|<f> (63) 

with T ((a^a)^) = 7/(0^0)^ . Or, way around, at ini- 
tial times the rate does not depend on the mean photon 



However both equations, (60) and (62) have the vacuum 
state as the asymptotic stationary state. Subtraction of 
(61) from (65) leads to 



, + \(E) I t \(<») -•viV^ (7*)' 

Wa)\ - (a?a)\ ' = e ' ^ — - 

1=1 



1 00 



Y np n+ i > 0, (66) 



77=1 



which shows that by using the nonlinear operators 
E_,E + , the calculated mean photon number of the re- 
maining photons in the cavity is always higher than when 
a,a^ operators are used. Because, with _E._E + -operators, 
the photons are always subtracted sequentially, one-by- 
one from any field state that does not contain the vacuum 
as one of its components, while by using a, a* operators, 
and depending on the field state, there exists the pos- 
sibility that the photon-leak occurs through the escape 
of more than one photon at a time. Thus, mean photon 
number reduction rate for i?-operators does not follow an 
exponential law, it lasts a longer time to have the cavity 
reaching a given photon mean photon number than in 
the amplitude damping model. We will touch this point 
again in the next section for a few specific examples of 
field states. 
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B. Solutions of the master equation 

The operator master equation (58) results in the fol- 
lowing infinite set of coupled equations for the diagonal 
elements p n = (n\p\n) of the statistical operator in the 
Fock basis: 



Pn = Pn+l - (1 - Sno) Pn, 



(67) 



where dot means the derivative with respect to the 'slow 
time' t = "ft. Making the Laplace transformation 



Pn(«) =£{Pn{t)} 



e ST p n {t)dT 



and remembering that 

£{i>n{T)} = sp n (s) - p n (0), 

we obtain from (67) the following set of algebraic equa- 
tions: 

Pn+l(«) = ( S + l)fn(«) -Pn(0), « > 1, (68) 

Pi(s) = sp o (s)-p o (0). (69) 
Their consequence is the relation 



p k (s) = s(s + l) k -%(s) (s + If" Vo(0) 



k-l. 



(,s + i) fe -V(o) 



■Pk-i(O). (70) 



Suppose that only the states with n < k were excited ini- 
tially. Then ph+i(r) = and consequently Pk+i( s ) = 0- 
In this case the consequence of (70) is 



Po(0) , Pi(0) , pa(0) 
Pol s J — r 



+ 



s(s + 1) s(s + l) 2 
Pfc(0) 



s(s + l) fe ' 
Taking into account the relations 

£{exp(— at)} = (a + s)~ , 
£{f™e- 4 } =m!(a + l)- m - 1 , 



(71) 



and the expansion 



1 



1 



1 



1 



1 



s(s + l) k s (s + 1) (s + 1) 2 (s + l) fe ' 

we can find the inverse Laplace transform of (71): 

Po(r) =Po(0)+pi(0) [l-e- T ] +p 2 (0) [l-e- T (l + r)] 
+ -"+Pfc(0)* fc _i(T), (72) 

where the function $fc(ar) was defined in (53). Having 
arrived at the expression (72), one can verify that it 
does not depend on the initial auxiliary assumptions that 
Pn(0) = for n > k, but it holds for any initial distri- 
bution. Knowing po(r) one can find all other function 



Pti(t) from equations (67). Finally, one arrives at the 
expression which coincides exactly with (56): 



Pm(r) = e Pm+k(0), m>l. 



(73) 



fe=0 



The reduced generating function of the diagonal ma- 
trix elements depends on time as follows, 



G(z;T) = X>"p„ = e -*X>(0)£ 



z n T k—n 



■i (k-n)\ 

n=l fe=l n=l v y 

Thus the mean number of photons evolves as 
dG 



n(r) 



dz 



' nT k ~ n 



z=l 



For the initial fc-photon Fock state \k) we obtain 



P 



(Fock) 



(r)=e" T 



-k—m 



{k-m)V 



1 < m < fc, 



T^ ocfe )(r) =e- T £] 



nr 



k—n 



^ (k-n)\ 

n=l x ' 



(74) 



(75) 



(76) 



(77) 



C. Special cases 

The series in the right-hand side of equation (73) can 
be calculated analytically for the initial negative binomial 
distribution (the corresponding pure negative binomial 
states were introduced independently in [56-58]) 



p(negbin) ^ _ T(/i + n) ^ U q 

T(fj,)n\ (n + M) 



7T. M>0, (78) 



where n is the initial average number of photons. The 
result is expressed in terms of the confluent hypergeo- 
metric function: 

P„(t) =p„(0)e- T $ (n + n;n+l;^-) (79) 
= p n (0) cxp (--^-)<S>(l-^n+l ;= ^-). (80) 



n + [i 



n + £t 



For integral values fj, = 1,2,... formula (80) can be writ- 
ten in terms of the associated Laguerre polynomials: 



Pn(r) = ,_ , " u+n cxp -— — - 



n +p, 



The special case of p, = 1 corresponds to the initial ther- 
mal distribution (coherent phase states in the case of 
pure quantum states) 



p„(0) =</(! + no) 



— \n+l 



(81) 
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In this case the photon number distribution preserves its 
form: 



P { t\r) =p„(0)exp[-r/(l + no) 



(82) 



and the mean number of photons decreases with time 
exponentially, although the rate of decrease diminishes 
with increase of the initial mean number: 



nW(T)=n exp [-r/(l + n )] 



(83) 



A formal substitution /j, — > — M with an integer M 
transforms the negative binomial distribution (78) to the 
binomial distribution (pure binomial states were consid- 
ered in [59] and rediscovered in [60,61]) 



P 



( '"")(0)= Mm o(M-n ) 



M-n 



n < M. (84) 



(M~n)\n\M M 
Then equation (79) is transformed to the formula 

M-n 



Pn(r) = ( - 



^ _ no 

M 



J M-n 



n T 
n -M 



In the case of the initial Poissonian distribution (co- 
herent pure states), 

n n 

Pn(0) = -5-exp(-Hb), (85) 
n\ 

the series (73) is reduced to the modified Bessel function: 



(Pots) 



(r) 



(no It) 



(86) 



OO 



n(Pois) {r) = e -n -r J- n (n /r) n/2 I n (2^) . (87) 

71=1 



The expressions (86)-(87) can be simplified in the asymp- 
totical case nor ^> 1, when the modified Bessel functions 
can be replaced by exponentials. Actually, one needs an 
additional condition t > n to ensure that the simpli- 
fied probabilities result in a convergent series whose value 
does not exceed 1. Thus for r 3> n + n^ 1 we can write 



P { n Pms) (r) 



(n /r) 



i/2 



■ exp 



(V%~Vt) 2 , (88) 



n {Pois) 



(r) 



tV4 



t 3 / 4 exp 



- (vW - Vt)' 



(89) 



In figure 4 we compare n(r) /n(0) for Fock, thermal and 
coherent states from (77), (83) and (87), respectively, 
with the corresponding result from amplitude damp- 
ing model (61). While the amplitude damping model 
presents an exponential decay independent of the field 
amplitude the exponential phase damping does not, the 
field taking longer to relax in this situation. Only for 



k = 1 and Fock states the two processes coincide. On 
comparing the curves for Fock, thermal and coherent 
states we observe that the decay rate shows a strongly 
dependence on the field state statistics. The more super- 
Poissonian is the field the longer will it take to decay 
in the model based on the exponential phase operators. 
It is interesting to note this dependence of the relax- 
ation process to the field statistics. In discussions of the 
amplitude relaxation process (e.g. [62,63]), this basis de- 
pendence is only attributed to coherence properties of 
quantum fields. In this aspect the coherent state is said 
to be selected from all the states as the more robust to 
decoherence from amplitude damping model [63]. We 
leave the discussion of coherence properties of the expo- 
nential phase model to a future publication, but we can 
anticipate from figure 4 that coherent states are not the 
more robust to dissipation in this model. 



V. CONCLUSION 

Summarizing, in this paper we proposed modifications 
in the SD photocount theory in order satisfy all the pre- 
cepts, as proposed by Srinivas and Davies for a consis- 
tent theory. Our central assumption was the choice of 
the exponential phase operators E_ and E+ as real 'an- 
nihilation' and 'creation' operators in the photocounting 
process, instead of a and . The introduction of those 
operators in the continuous photocount theory, besides 
eliminating inconsistencies, leads to new interesting re- 
sults related to the counting statistics. A remarkable 
result, which is responsible for all the physical consis- 
tency of the model, is that in this new form an infinitesi- 
mal photocount operation J E really takes out one photon 
from the field, if the vacuum state is not present. Con- 
sequently the photocounting probability distribution for 
a Fock field state is Poissonian, evidencing again the di- 
rect correspondence of number of counted photons and 
number of photons taken from the field. 

We also have investigated the evolution of the field 
state when photons are counted, but with no readout, 
leading to the pre-selected state. The mean photon num- 
ber change shows now (in contrast to the exponential 
law obtained for the amplitude damping model) a non- 
exponential law, which only depends on the condition 
that photons are present in the field, independently of 
their mean number. 

An advantage of the proposed model is its mathemati- 
cal consistence. Since many of its predictions, especially 
those related to multiphoton events, are significantly dif- 
ferent from the predictions of the SD theory, it can be 
verified experimentally. One of the first questions which 
could be answered is: whether the decrease of number of 
photons in the cavity due to their continuous counting 
always obeys the exponential law (61) (i.e., the rate of 
change is proportional to the instantaneous mean num- 
ber of photons), or noncxponential dependences can be 
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also observed (for example, in the case of detectors with 
large dead times)? 

We leave for a future work a detailed study of pho- 
tocounting processes in the presence of other incoher- 
ent (dissipative) processes and discussion about coher- 
ence properties of the field under the exponential phase 
damping model. Also, an open question is how extremely 
nonlinear (with respect to a and a?) one-count operator 
(25) could be derived from some microscopical model, 
starting from fundamental interactions which are linear 
with respect to operators a and a) (as soon as the model 
considered in our paper is pure phenomenological, as well 
as the SD model). 

We would like to emphasize that many preceding stud- 
ies, such as [11-16], adopted the SD theory of photode- 
tection. We believe that applications of our consistent 
model of photodetection to these problems may, indeed, 
bring new and important results for both the quantum 
measurement theory and experiment. 
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Figure captions 



Figure 1. Conditional photocount probability distribution for the initial number state with m — 5 photons. Solid 
lines are for the present model while dashed ones are for the original SD theory. Numbers above the curves correspond 
to the fc-event. 



Figure 2. Same as figure 1 for a coherent state with average \a\ 2 = 5 photons. 



Figure 3. Same as figure 1 for a thermal state with average n = 5 photons. 



Figure 4. Normalized mean number of photons n(t) = n(jt)/n(0) in the cavity under continuous measurement, 
for mean initial number of photons n(0) — 1,5 and 10. Different line styles represent different states for each mean 
initial photon number. Three lower curves correspond, in the order from bottom to top, to the Fock state (solid line), 
coherent state (dashed), and thermal state (dotted), for n(0) = 1. Three middle curves are related to the case of 
n(0) = 5 in the following order (from bottom to top): Fock state (dash-dotted line), coherent state (dash-dot-dotted), 
thermal state (short dashed). Three upper curves are related to the case of n(0) = 10 in the same order as before: 
Fock state (short dotted line), coherent state (short dash-dotted), thermal state (solid). The lowest solid line (the 
Fock state with n(0) = 1) coincides with the exponential decay for the amplitude damping model. 
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Fig. 2 - M.C. de Oliveira et al. 
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Fig. 4 - M.C. de Oliveira et al. 




